THE NEWTON WONDER IN MECHANICS 



By D. Lynden-Bell 
Institute of Astronomy, Cambridge 

Application of Newton's ideas from Principia gives many new results 
in mechanics. Here we explore the question "What form of extra force 
will maintain the magnitude of a vector constant of the motion while 



(N 



changing its direction 



Introduction 

Shortly after his death I reviewed Chandrasekhar's last book 1 Newton's Principia 
for the Common Reader. This fine book introduced me to the whole sweep of Newton's 
ideas in mechanics and sent me back to reading Principia itself in translation 3 . New- 
ton's insight stimulated many new ideas in me of which only the following have been 
§ ! investigated so far 

^ : 

Revolving orbits 

O , 

In his theorum on revolving orbits Newton points out that in a central orbit r 2 d(f)/dt 
h then r 2 d(acf))/dt = ah. So an angular motion with 0* = a<p replacing will still 



have constant angular momentum. He then asks what extra central force is needed to 
leave the r motion unchanged and deduces that it must balance the extra centrifugal 
force (a 2 — l)h 2 r~ 3 . He then deduces that the apses of such orbits rotate at the rate 
(a — 1)2tt/P faster. Here P was the mean period in of the original orbit. The new 
orbit when viewed from axes revolving at the rate (a — l)d<p/dt is exactly the same as 
the old one viewed from fixed axes. But (a — l)d<fi/dt is not a steady rotation rate. In 
reference 4 we explore the weird distorted shapes taken by such orbits when viewed from 
axes rotating at the mean rate (a — 1) < d<p / dt > where < d<p/dt >= 2% /P. 
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Exact general N-body solutions 

Few know that Newton solved completely an N body problem for any N for any 
initial conditions. The force law for which he did this was 

F 7 j = km I mj(r J - 17). 

The total force on body / is then directed at the barycentre at r 

F I = J2 F/j = kmjM(r - n) 

where M is the total mass. Thus each body moves as though attracted to the centre of 
mass which itself travels uniformly. This force law can be found by differentiation of the 
total potential energy V = \kM 2 r 2 where 

r2 = J2J2 m i m J M ~ 2 ( r J ~ T if- 

i< J 

Ruth Lynden-Bell and I have shown 5 that this result can be generalised to solve the 
non-harmonic case in which V is replaced by any V(r). Some results even extend to V 
of the form 

where V is an arbitrary function of r and V 2 is any function of its many scale-free 
arguments. After some struggles with Einstein-Bose & Fermi-Dirac statistics we solved 
the V(r) problem completely in quantum mechanics 6 . We have further work in progress 
on the statistical mechanics and mechanics of these extraordinary N-body problems 
which can be applied to Bose-Einstein condensates. 

Vector constants of the motion 

We all know that in the presence of a central force the radius vector from the centre 
to a body sweeps out equal areas in equal times in a plane. However, although Newton 
proves it, we are not taught the converse that if from some fixed point S, the radius 
vector to a body sweeps out equal areas in equal times in a plane then the force on 
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that body is central to S. This is Newton's 2 proposition 2, but he follows it with a 
truly remarkable scholium or rider that no-one seems to have understood. If we have 
a body subjected to a central force but beside that there is another force which is at 
each moment perpendicular to the current plane of motion, then the radius vector to 
the body still sweeps out equal areas in equal times, notwithstanding its non-coplanar 
motion. 

To prove this consider the moment of this extra force about the centre. Since the 
extra force is parallel to the angular momentum the moment is perpendicular to it. Thus 
the rate of change of the angular momentum is perpendicular to it, so the magnitude 
or 'length' of the angular momentum vector remains constant. In the surface made up 
of infinitessimal areas perpendicular to the angular momentum at each instant the rate 
of sweeping of area ±r 2 d<f)/dt = h/2 is therefore constant QED. This is all so obvious to 
Newton that he does not stoop to prove it - he simply states it! 

Analysing what Newton has done here, he takes a vector constant of the motion 
and asks under what circumstances the magnitude of that vector can be conserved 
without requiring each component to be separately conserved. But there are three vector 
conserved quantities in classical mechanics. They are, besides the angular momentum, 
the linear momentum p and r — vt = r , the position of the barycentre at t — 0. Can we 
preserve |p| without keeping p constant? Yes it requires a force dp/dt = mv x B that 
is always perpendicular to p. Thus we are led to forces of gravo- magnetic type of which 
Coriolis force is a special example. What forces are needed to preserve the magnitude 
of r — vt? Evidently (r — vt) • J^(r — vt) = so — vt must be of the form (r — vt) x B. 

Thus the force per unit mass F must be of the form 



where B may be any function of position, time and velocity. Perhaps the simplest 
example of such a force occurs when a magnetic field is constant in space but varies 
its strength with time. The associated electric field is E = ±r x dH/dt so the total 
electromagnetic force on unit charge with unit mass will be of the form (1) provided 



F = vxB-rx B/t 
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B = bt 2 with b a constant. We orient the z axis along it so b = (0, 0, b). Using the 
fact that (r — vt) 2 = const = z 2 + Rq the equations of motion may be solved giving 

x + iy = t(u + we tb / t ) ; x — v x t + i(y — v y t) = ibwe lb ^ 

V J ) (2) 

z = v z t + z 

where u and w are complex constants and v z & z are real ones. \w\ = Ro/b. Thus New- 
ton's method of generalising constants of the motion works beautifully for this problem 
too. 

Special relativity 

Can the idea be extended to relativity? Here the ten classical integrals combine into 
the energy-momentum 4-vector, P and the anti-symmetric tensor 

K-ij tijkix P . 

The invariant scalars are 

P 2 , KijK ij = 4K 2 and e ijkl K i:j K kl . 

Unfortunately P 2 = m^c 2 which is always constant and the epsilon KK invariant is zero 
identically, so only K 2 is suitable for playing Newton's game. Differentiating with respect 
to comoving time r and remembering that mrfx/rfr = P, we find that K 2 is constant 
when the four-force F = dP/dr is perpendicular to x = (ct, x, y, z). However, since 
P 2 is constant the four-force is always perpendicular to P. If N is any vector along a 
third independent direction perpendicular to the four-force, that force takes the form 
F% — £ijkix J P h N l so this is the general form of force when K 2 = (r x p) 2 — E 2 c~ 2 (r — vt) 2 
is conserved. Here E is the energy. A pretty example is given by taking N along the 
time direction; then reverting to 3- vectors and writing N — |N| 

dp/dr = JVrxp. (3) 

Under such forces |p| is clearly constant so |v| = v is constant too; furthermore writing 
L = r x p we see that dh/dr = Nr x L so |L| is constant. From our K 2 invariant this 
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of course implies |r — vt\ is constant, p and L are perpendicular and obey the same 
equation (3) which corresponds to angular rotation at the variable rate Nr^Jl — v 2 /c 2 . 
A third vector perpendicular to the other two is given by R = p x L/p 2 which is the 
component of r transverse to p. It is readily shown that it too obeys dH/dr = iVrxR, so 
|R| is constant and R, p, L form a triply orthogonal triad of rotating vectors of constant 
length, v is along p and is also fixed in that rotating frame. Now 

d(r ■ v)/dt = v 2 = const, 

so 

r = R + v(t + a) (4) 

since the components of r transverse to v make up R. The above equation gives the 
misleading impression that the particle is travelling uniformly in a straight line. That is 
indeed true relative to the rotating axes but in absolute axes R and v are rotating 
and at a time dependent rate. Notice that X = r — vt is constant in the rotating axes 
so X moves on a sphere in fixed axes i.e., |X| = const. Viewed from the rotating axes 
in which R, p and L are fixed constant vectors the whole motion is simple and in those 
axes even their rotation rate takes the simple form, 

= (1 - v 2 /c 2 )^Nr = (1 - v 2 /c 2 )^N [R + v(t + a)} (5) 

although it varies in magnitude and direction. Of course from the viewpoint of a co- 
rotating observer the universe rotates at rate —ft. The true complication of the motion 
is only comprehended when one imagines what is seen by an inertial observer to whom 
the inclined offset straight line, down which the motion occurs uniformly, whirls around 
a changing axis at a non-uniform rate. So far the reader may have assumed that N was 
a constant but nowhere have we assumed this and indeed the force law can be such that 
N is any scalar function of r, v and t. 

In parting we notice that this is but a very special case of the force laws that preserve 
K 2 . We chose it because the others necessarily involve t explicitly but it is peculiar in 
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that N's direction is fixed (along the time direction) and not varying as it would in the 
general case. (3) is in fact a special case of (1) with B radial. 



The eccentricity vector and magnetic monopoles 
To return to where Newton started, his equation of motion for a particle moving 
under a central force — V'r with an extra force iVr x v reads for unit mass 

dv/dt = -V'r + Nr x v (6) 

with the central force present can we still find some rotating axes in which the motion is 
especially simple? Letting dots denote rates of change with respect to the moving axes: 
dv/dt = v + Q x v. Evidently we should again take fi = iVr. Then dr/dt = r = v so 
even when N varies 

r = -V'r (7) 

so if we take the rotating axes the whole dynamics is the same as if the N force were 
absent and we had fixed axes. Thus in the rotating axes we have r x r = h = 
constant relative to the rotating axes and when V = —GM/r we write r = r/r and 
obtain 

h x r = —GM (r x *) x r = —GMi 



so 

hxr= -GM(r + e) 
where e is a vector constant of integration, e is constant in the rotating axes not 
the fixed ones. Those with no history call it the Runge Lenz Vectorf] but to me it 
is Hamilton's eccentricity vector pointing to pericentre with magnitude equal to the 
eccentricity. Evidently 

e = r x h/(GM) - r = v x h/(GM) - r (9) 



Hamilton 7 gave the first coordinate free derivation in 1845 by writing the equations in quaternion 
form. Bernoulli 8 , Laplace 9 and possibly even Newton got the result earlier. Atomic Physicists learned 
the result from Runge & Lenz much later. See Goldstein 10 
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again e is constant in our rotating axes, h is not constant in the fixed axes but obeys 
dh/dt = Nrxh. 

We are now in a position to ask a question that crossed my mind but which I would 
not have remembered to investigate had not S Aarseth reminded me. What is the most 
general force law that leaves |e| invariant while letting e change? 

Evidently the example given above has this property since e precesses with the 
axes at the rate ft = Nr which may well be variable in magnitude as well as direction; 
however that example does not give the general case. We now revert to fixed inertial 
axes. Again h = r x v whether that quantity is constant or not and e is defined by (9) 
with v written for r. 

Differentiating GMe we have from (9) 

„. r de dv ( dv\ „, r dx 

GM— = — xh + vx r x — - GM— . 
dt dt \ dt J dt 

The extra force per unit mass that causes the change in e must be 
F = g + GMr- 2 r, so 

de 

GM— = Fxh + vx(rxF) (10) 

(JjL 

We want this extra force to leave e 2 unchanged so there must be no component of de/dt 
along e. This condition may be put in the form 

¥■£ = (11) 

where 

£ = hxe + (exv)xr . (12) 
Thus the general F that leaves e 2 invariant is of the form 

F = kx£ (13) 

and it will cause a change of e given by 

GM de/dt = (k x t) x h + v x (r x (k x £)) (14) 
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since the last term of £ itself involves two cross products we see that the final term above 
involves five consecutive cross products and it takes some work to sort them out and to 
re-express the RHS as a cross product with e so that 

de/dt = uj x e (15) 

We write 

C = (k • r)e + (k • e)r 
and notice it is in the plane of motion and find 

GMuj = 2(k • h)h + (v • r) (C x v) x e/e 2 (16) 

Thus our result is that the general force is of the form (13) in which k is any vector 
(which may vary with time) while it leaves e 2 invariant. The angular velocity of e is 
given by (16). The first term causes e to swing in the plane of motion while the seconds 
swings it down through that plane (C x v is parallel to h). A useful final check on 
the complicated derivation is that if any multiple of £ is added to k it must make no 
difference. Indeed C x v changes by (£ ■ r e + £ ■ e r) x v which is identically zero by 
(12). 

The general case discussed above requires complicated forces of the form given by 
(12) and (13) but a great simplification occurs if we demand that these forces do no work. 
The condition F • v = leads via (12) & (13) to k • h = and since we can always take 
k perpendicular to £ we set k = N£ x h/£ 2 . Then the extra force F reverts to Newton's 
simple form, Nr x v discussed earlier (6). We may consider this as a magnetic force by 
endowing our particle with unit charge and writing B = — cNr . It is then natural to 
require that divB = so r~ 2 d/dr(Nr 3 ) = and therefore cNr 3 = —Q = constant. The 
resulting field is B = Qr/r 2 which is clearly the field of a magnetic monopole of strength 
Q. As Poincare 11 first found, motion in such fields has a singular beauty. For a particle 
of mass m and charge q orbiting a monopole of charge — Zq and monopole strength Q 

mr = —(rr~ 2 + rjr x vr~ 2 
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where ( = Zq 2 and rj = Qq/c. Cross multiplying by r, dL/dt = — r/h x r/r 2 = 
— 77 (r x v) x r/r = —77 dr/dt where L = mh. It follows that L 2 is constant and that 

L — r/r = J = const. 

From which it follows that L 2 = J • L and rj = — J • f so L and r precess around J on 
coaxial cones. The orbit is an ellipse that precesses around the r cone because of the 
deficit angle that is missing when the cone is flattened, but I gave the details earlier 
and solved the quantum mechanics of monopolar hydrogen. The results here can be 
extended to give the interesting motions of a small spinning sphere both classically and 
in quantum mechanics 12 . Even the Dirac equation has a pretty solution in the field of a 
charged monopole. In general relativity the equivalent of a monopole is a gravomagnetic 
monopole and the required generalisation of Schwarzschild's metric is called, NUT space, 
it has many weird properties. Although the space has the same properties when viewed 
from any direction its metric is not spherically symmetric and can not be transformed 
into spherically symmetric form 12 . Thus following a lead from Newton one ends up 
challenging the concept of spherical symmetry in General Relativity! 

Nouri-Zonoz and I have explored the gravitational lensing of gravomagnetic monopoles 
13 and he has found the cylindrically symmetrical space corresponding to a line gravo- 
magnetic monopole 14 - the metric can not be put into a single-valued cylindrically 
symmetric form! 

What we have found so far comes from a tiny fraction of Principia; we hope that 
it will encourage others to read Chandrasekhar's last book and Newton's great one! 
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